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ABSTRACT 

A  situation  is  considered  in  which  a  grenade  launcher  fires  at  a  target, 
observes  the  projectile  impact,  adjusts  the  point  of  aim,  fires  a  second 
round,  adjusts,  fires,  and  so  on.   A  model  is  developed,  based  on  the 
assumptions  of  normal  ballistics,  perfect  observation  of  impacts,  adjust- 
ments without  error,  and  a  unimodal  target  destruction  function.   The 
problem  is  to  determine  optimal  adjustments,  in  order  to  maximize  the 
probability  of  target  destruction  within  a  given  number  of  rounds.   It 
is  shown  that  seemingly  different  adjustment  procedures  are  equivalent, 
if  viewed  in  appropriate  coordinate  systems.  Previous  results  concerning 
sequential  adjustments  which  are  constrained  to  be  linear  functions  of 
observed  impact  points  are  extended  to  the  class  of  translation  invariant 
procedures.  Properties  of  the  optimal  sequential  adjustment  procedure, 
including  some  related  to  stochastic  approximation,  are  reviewed.   The 
effects  of  errors  in  judging  impact  positions  are  discussed. 

Prepared  by: 


1.   Introduction. 

Consider  the  problem  of  adjusting  grenade  launcher  fire  on  a  target, 
using  a  "shoot-adjust-shoot"  strategy.  After  each  round  is  fired,  the 
grenadier  observes  the  projectile  impact,  and  estimates  its  distance  from 
the  target.  He  then  attempts  to  adjust  his  point  of  aim,  so  as  to  place 
the  mean  of  the  ballistic  distribution  of  impact  points  directly  over  the 
target.   The  problem  we  address  is  that  of  determining  an  optimal  adjust- 
ment scheme. 

In  what  follows,  we  assume  the  ballistic  distribution  of  impacts 
of  rounds  fired  at  a  given  aimpoint  is  bivariate  normal.  We  further 
assume  that  the  target  is  a  point  target,  and  that  the  target  destruction 
function  (or  lethality  function)   d  is  the  same  for  all  rounds  fired. 
The  value  of  d(x,y)   is  defined  to  be  the  probability  that  the  target 
is  destroyed  by  a  round  impacting  at  the  point   (x,y)   relative  to  the 
target.  The  difficulty  of  finding  optimal  adjustment  procedures  hinges 
initially  upon  how  errors  are  made  by  the  grenadier  in  estimating  the 
position  of  each  impact  point,  relative  to  the  target.   For  the  most  part, 
we  confine  ourselves  to  the  simplest  case  in  which  the  grenadier  is 
"perfect" — he  delivers  adjustments  which  are  exactly  the  desired  propor- 
tions of  the  observed  miss  distances.   In  the  last  section,  some  remarks 
are  made  concerning  the  probable  effects  of  errors  in  this  process. 

The  problem  of  how  to  optimally  adjust  fire  on  a  stationary  point 
target  has  been  of  considerable  interest.  Grubbs  [8]  gives  an  adjustment 
scheme  which  is  "best"  (i.e.  minimum  variance)  and  unbiased  among  linear 
functions  of  the  impact  points.   (In  the  literature  of  statistical 
estimation,  such  a  scheme  is  called  "BLUE",  for  "best  l_inear  unbiased 


estimator .")  Nadler  and  Eilbott  [10]  derive  the  BLUE  adjustment  scheme 
and  show  it  maximizes,  among  linear  functions  of  the  impact  points,  the 
probability  of  destroying  the  target  with  each  round,  for  a  fairly  general 
class  of  destruction  functions.  This  procedure  is  a  special  case  of 
stochastic  approximation  schemes,  about  which  a  great  deal,  mostly  of  an 
asymptotic  nature,  has  been  written  [5],  [11],  [12].  Other  stochastic 
approximation  approaches  to  the  adjustment  of  fire  have  included  cases 
in  which  the  observer  cannot  precisely  locate  the  impact  points  (see, 
for  example,  Barr  and  Piper  [3]). 

The  main  purpose  of  this  paper  is  to  demonstrate  that,  with  appro- 
priate coordinate  systems,  the  procedures  of  Grubbs  and  of  Nadler  and 
Eilbott  may  be  obtained  from  well  known  statistical  theory  concerning 
BLUE  estimators.  We  also  show  (in  Section  2)  that  these  procedures  are 
indeed  stochastic  approximation  procedures,  and  (in  Section  3)  that  they 
possess  an  optimality  property  stronger  than  BLUE.   The  present  approach 
seems  notationally  and  conceptually  simpler  than  those  used  previously. 

It  is  notationally  convenient  to  restrict  our  consideration 
to  a  single  dimension,  say  range,  in  what  follows.  Under  the  assumptions 
we  make  below  (which  are  somewhat  weaker  than  those  of  Nadler  and  Eilbott) , 
the  restriction  to  one  dimension  causes  no  loss  of  generality:   identical 
optimal  one-dimensional  procedures  carried  out  independently  in  each 
dimension  gives  an  optimal  two-dimensional  procedure.   (See  Section  5  below.) 

2.   The  Model. 

Assumptions  which  permit  analysis  of  a  two-dimensional  problem  by 
separate  analyses  with  each  dimension  are  essentially  those  required  for 
Anderson's  [2]  inequality.  The  following  are  sufficient: 


i)   statistical  independence  of  the  range  and  deviation  components 
of  the  impact  points,  each  component  having  a  symmetric  uni- 
modal  continuous  distribution, 

ii)   no  grenadier  errors  are  made,  and 

iii)   the  destruction  function  d  is  symmetric  about  the  impact  point 
(x,y)   and  {(x,y):  d(x,y)  k  a}  is  convex  for  all  a  ^  0. 

In  particular,  a  circular  normal  ballistic  distribution  with  a 
destruction  function  of  the  "circular  normal"  or  "circular  cookie-cutter" 
types  will  satisfy  condition  (iii) .   It  is  well  known  that  in  actual 
gunnery  problems,  the  ballistic  distribution  is  such  that  the  variance 
in  range  is  not  the  same  as  that  in  deviation.   In  most  gunnery  problems, 
these  variances  depend  on  many  factors,  but  principally  upon  the  grenade 
launcher-target  range.   The  assumption  that  the  impact  points  fall  in  a 
circular  normal  pattern  about  the  aimpoint  is  equivalent  to  agreeing  to 
measure  range  and  deviation  in  the  corresponding  standard  deviation  units, 
and  in  addition,  assuming  independence  between  the  range  and  deviation 
components  of  the  impact  points. 

Suppose,  then,  that  XJ  ,X*  ,xl: , . . .   are  independent  normal  random 
variables  with  unknown  mean  u  and  known  variance  a2   (we  may  take 
a^  =  1  without  loss  of  generality  ) .   These  random  variables  represent 
the  sequence  of  impact  point  ranges  that  would  be  obtained  with  a  sequence 
of  rounds  fired  without  adjustment  (i.e.  with  a  fixed  aim  point).   The 
problem  is  to  determine  a  sequence  of  adjustments,   £-t.t€9i*<*i  where  £ 
is  the  adjustment  made  between  the  n —  and  n  +  1 —  rounds,  so  that 
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As  mentioned  above,  assuming  a*-   «  1  is  equivalent  to  measuring  distances 

in  a -units.   In  theory,  a2  varies  with  the  range  being  fired,  and  hence 

varies  as  adjustments  in  range  are  made.  However,  in  practice,  this 

variation  in  a2  is  sufficiently  small  within  a  given  firing  mission  that 

a2  may  be  assumed  to  be  constant. 


the  probability  of  target  destruction  with  each  round  is  maximized.   (As 
we  remark  below,  this  will  maximize  the  probability  of  target  destruction 
within  any  fixed  number  of  rounds.) 

In  general,  the  n —  adjustment  is  a  function  of  the  first  n 
observed  impact  ranges,  X..,X9,...,X  .   Each  adjustment  £  has  the 
effect  of  translating  the  coordinate  system  -£  units.   For  example,  if 
an  adjustment  £  =  -50   (corresponding  to  the  intention  of  dropping  50 
meters  in  range)  is  made,  the  net  effect  is  as  if  the  coordinate  system 
(i.e.  target)  were  transformed  50  units  away  from  the  grenade  launcher. 
We  must  therefore  keep  track  of  the  adjustments  made  up  to  the  n — 
round,  so  as  to  be  able  to  express  the  n —  impact  point  in  coordinates 
which  are  meaningfully  related  to  the  previous  impact  points.  Accordingly, 
let  X.   denote  the  j —  impact  range,  expressed  in  the  r —  coordinate 
system.  Upon  firing  the  first  round,  the  grenadier  observes  x| ,   and 
makes  adjustment  £..  .   The  first  round  can  thus  be  given  as  X^  =  X*  +  £-, . 
Similarly,  the  second  round  X^,   fired  after  the  adjustment  £-.  ,   can  be 
given  in  the  initial  coordinate  system  as  X*  ■  X^  -  £_ ,   and  so  on  for 
subsequent  rounds,  adjustments  and  coordinate  systems. 

Remark.   It  should  be  noted  that  the  procedure  we  will  determine 
does  not  require  the  grenadier  to  actually  estimate  the  miss  distances  X.. 
Rather,  it  will  turn  out  that  he  should  attempt  to  adjust  by  a  fraction 
of  the  observed  miss  distance.   Thus,  the  grenadier  needs  only  to  be  able 
to  duplicate  some  fraction  of  an  observed  distance  in  his  aim  adjustment. 
The  present  discussion  involving  estimates  X   of  actual  distance  values 
are  a  theoretical  convenience  only. 


If  all  observed  impact  ranges  are  transformed  to  the  initial 
coordinate  system,  and  if  each  adjustment  is  some  linear  function  of  the 

data  available  for  that  adjustment  (i.e.   £  (X*  ,X2 X  )),  then  the 

£'s  are  also  linear  functions  of  the  "unadjusted"  impact  ranges 
xj ,X* . . . ,XX .  The  problem  then  becomes  that  of  finding  the  optimal 
linear  function  of  a  random  sample  of  identical,  independent  normal 

(y,l)   random  variables.  Naturally,  one  expects  the  sample  mean 

1  n 
1  =  —  Y  x} ,   to  provide  an  optimal  adjustment  -C1 ,   from  the  initial 


n 
X 

aimpoint  u,  where  we  write  "51,"  rather  than  "£  "  to  indicate  we 


mean  adjustment  from  the  initial  aimpoint.   In  what  follows,  we  call  this 

the  "X-bar  procedure." 

It  is  easy  to  establish  that,  in  fact,   C1  =  -X1   is  the  optimal 

linear  adjustment.   It  is  well  known  that  X1   is  a  BLUE  estimator  for 
J  n 

u  under  the  present  conditions.   With  the  unimodal,  symmetric  target 
destruction  function  considered  here  it  is  clear  that  one  may  restrict 
his  attention  to  linear  adjustments  which  are  unbiased  estimators  of  -y, 
since  otherwise  the  probability  of  target  destruction  on  the  next  round 

could  be  made  as  high  or  higher  by  subtracting  the  bias.   But,  among  the 

st 
normally  distributed  aimpoints  for  the  n  +  1 —  round,  which  arise  after 

applying  unbiased  adjustments  51,   the  probability  of  target  destruction 

is  maximized  by  minimizing  the  variance  of  the  aimpoints  about  the  target. 

That  is  precisely  what  the  BLUE  estimator  of  -u  achieves,  and  it  is  the 

optimal  scheme  developed  by  the  aforementioned  authors. 


For  example,  Larson  [9],  gives  the  following  [Theorem  7.3.1]:  X  is  a 
random  variable  with  mean  u  and  variance  a2.   If  X-,X2»...,Xn  is  a 

random  sample  of  X,  X  is  the  best  linear  unbiased  estimator  of  u. 


In  order  to  see  that  the  X-bar  procedure  {E1}   is  a  stochastic 

n 

approximation  procedure,  we  shall  express  it  in  terms  of  the  sequence  of 
current  coordinate  systems,   {£.}.  Note  that 


n   xi   V 


After  making  this  adjustment  and  firing  the  second  round,  resulting  in 
X^,  we  have,  in  the  second  coordinate  system,  two  observations:  X^  and 


X* .  Now 


r2  =  pi  _  pi 
t,2    i>2         ^ 


=  ~  2  (X2+Xi)  '  ?1 


±X* 
2  V 


In  general,  £  =  E1  =  51  . 
n    n    n-1 


-  -i  (Xl+(n-l)^  J  -  CI  , 
n   n      n-1     n-1 


=  -  1  xn. 
n  n 

Thus,  the  optimal  "X-bar"  procedure  is  identical  to  the  "  

stochastic  approximation"  [11]  procedure  which,  after  the  n —  round, 
adjusts  by  the  fraction  —  of  the  observed  miss  distance  of  the  n — 
round.  The  good  asymptotic  properties  of  such  stochastic  approximation 
procedures  are  not  surprising  when  viewed  in  this  light.   Similarly, 
the  distribution  of  the  n  +  1 —  impact  point  X  . ,   derived  in  [10] 
using  multivariate  analysis,  is  now  very  easily  seen  to  be  normal  with 


mean  zero  and  variance   (1  +  — )a2  ■  - .  That  the  mean  ia  zero  follows 

n       n 

from  the  fact  that  X   is  unbiased  for  u  -  E(X*   ).  The  variance  of 

n  n-t-i 

the  aimpoint  for  the  n  +  1—  round  is  VQC1)  =  -,   and  since  X^  is 
independent  of  X1 ,  V(X^)  =  VO^)  +  V(X*)  -it^ 

Remark.   It  is  easily  shown,  by  induction  for  example,  that  the 

n—  adjustment  in  terms  of  the  j —  coordinate  system  (data  and  shift) 

1  n   ' 

is  that  which  places  XJ  at  zero  (i.e.  on  the  target);   £  =  -  —  I     X. . 

n  n     n  ^_^  i 

The  earlier  observation  that  £n  =  Xn  also  follows  from  this,  since 

n     n  n 
n-1 

I  X»E0. 
i=l 

3.   Strong  Optimality  of  the  X-bar  Procedure. 

We  have  re-established,  by  elementary  methods,  that  the  X-bar 

procedure  has,  under  our  assumptions  of  independence,  normality  and 

regularity  of  the  target  destruction  function,   the  properties 

(a)  among  linear  adjustments,  it  maximizes  the  probability  of 
destroying  the  target  with  each  round  fired,  and 

(b)  it  provides  unbiased  aimpoints  (after  the  initial  round  is  fired) 
with  minimum  variance.   It  easily  follows  that  also  it 

(c)  maximizes  the  probability  of  destroying  the  target  within  the 
first  k  rounds;  k  =  1,2,...,  and 

(d)  minimizes  the  expected  number  of  rounds  required  to  destroy  the 
target.   The  latter  may  be  established  using  the  identity  [4] 

00 

E(N)  =  I     P[N  i  n], 
j-l 
where  N  is  the  (random)  number  of  rounds  required  to  destroy 
the  targeUEvans  [6]  gives  an  elegant  proof  of  the  intuitively 
obvious  fact  that  maximizing  the  probabilities  of  target  destruc- 
tion with  each  round  minimizes  the  expected  number  of  rounds 
required  to  destroy  the  target. 

Since  X  is  a  robust  estimator  of  y,  with  central  limit  theorem 
properties  when  the  X.'s  are  not  necessarily  independent  and  identical,  it 
follows  that  the  X-bar  adjustment  also,  in  loose  terms 
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(e)  is  robust  (i.e.  tends  to  operate  efficiently  even  when  the 
underlying  assumptions  are  not  met) ,  and 

(f)  tends  to  give  normally  distributed  aim  points. 

We  wish  to  show  next  that  the  X-bar  procedure  is  optimal  in  a 
class  of  procedures  much  wider  than  the  class  of  linear  procedures 
considered  heretofore.   The  approach  is  to  use  a  result  of  Wolfowitz  [3], 
concerning  minimax  estimators,  in  a  decision  theoretic  setting.   We  use 
the  notation  and  terminology  of  Ferguson  [7] .   In  the  preceding  section, 
it  was  established  that  the  problem  of  finding  an  optimal  aimpoint  after 
n  rounds  have  been  fired  is  equivalent  to  finding  an  optimal  estimator 
for  y,   based  on  a  random  sample  x|,X^,...,X^   from  a  N(y,l)   popula- 
tion; once  that  estimator  u  is  found,  the  optimal  adjustment  for  the 

next  round  is  -u,   from  the  initial  aimpoint.   Clearly,  this  problem  is 

* 
invariant  under  the  group  of  translations   (i.e.  roughly  speaking,  if  y 

is  optimal  for  y,   y  +  c  is  optimal  for  y  +  c) .   It  is,  therefore, 

extremely  reasonable  to  insist  that  the  solutions  to  the  problem  should 

be  invariant  (i.e.  roughly  speaking,  if  y(XJ,...,X1)   is  optimal  for  y, 

then  y(x}+c,X;!;+c, . . .  ,Xx+c)   should  be  optimal  for  y  +  c)  .  We  shall, 
L  Z  n 

therefore,  restrict  our  attention  to  estimators  which  are  invariant  under 
translations;  this  means  essentially  that  we  are  indifferent  as  to  where 
the  origin  of  our  coordinate  system  is  placed,  relative  to  the  target, 
before  firing  begins.   (It  is  intuitively  unacceptable  for  an  arbitrary 
assignment  of  this  sort  to  affect  the  procedure  to  be  followed  in  adjust- 
ing fire.) 


For  a  discussion  of  invariant  problems  and  solutions,  see  Ferguson  [7], 
Definitions  4.1.2  and  4.2.1.   It  should  be  noted  that,  while  functional 
equations  of  the  form  f(x+c)  =  f(x)  +  c  determine  linear  functions  in 
the  single  variable  case,  this  is  certainly  not  the  case  for  functions 
of  several  variables  (see  Acz£l  [1]). 


Theorem.   In  the  class  I   of  translation  in  variant  procedures,  the 
X-bar  procedure  is  uniformly  best,  in  the  sense  that,  for  any  t  >  0, 

max  P[|y(xJ,...,X1)  -  y|  <  t] 

ie         1 

=  pMx1  -  ul  <  t]. 

1  n    ' 

Our  proof  makes  use  of  results  of  Wolfowitz,  which  are  stated  in 
convenient  form  by  Ferguson  (Theorem  4.6.1)  as  follows: 

Lemma.   If  x] , . . . .X1   Is  a  sample  from  a  normal  distribution  with  mean  y 
In 

and  known  variance  [1,  say],  then  X1   is  a  best  invariant  estimate  of  y 

and  a  minimax  estimate  of  u,  provided  that  the  loss  function  L(y,y) 

is  a  nondecreasing  function  of   lu  -  u|   and  that  EtLCOjX1)^  =  0] 

n 

exists  and  is  finite. 

To  say  that  X   is  a  best  invariant  estimate  of  y  means  that 
n 

it  is  as  good  as  any  other  invariant  estimate:   for  every  y  and  any 

u  €  I, 

E  L(y,xi)  s:  E  L(y  ,y(Xx , . . .  ,X*)) , 
y     n    y      in 

where  "E  "  denotes  expectation  with  respect  to  the  joint  normal  distri- 
bution in  which  the  X^   have  mean  u.  To  say  that  X1   is  minimax 
estimate  of  y  means  that,  roughly  speaking,  it  minimizes  the  maximum 
loss  we  could  sustain,  regardless  of  "natures  choice"  of  y. 

Now  let 

/  0  if   |y  -  y|  <  t 
L(y,y)  =  < 

1  otherwise. 
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Clearly  L  satisfies  the  conditions  of  the  lemma,  since  for  any  fixed 

t  >  0, 

i)   L  is  nondecreasing  in   |y  -  y|,   and 

ii)   E.LCO.X1)  =  P[|xM  ^  t|v  =  0]  £  1. 
U     n       n 

By  the  lemma,  it  follows  that,  for  any  U, 

E^L(y,y)  =  P[|y  -  u|  it] 

is  maximized  over  I   by  taking  y  =  X*  ,  which  completes  the  proof. 

The  theorem  may  be  interpreted,  along  the  lines  of  the  argument 

in  Section  2,  as  guaranteeing  that  the  X-bar  procedure 

(g)   has  properties  (a)  through  (d)  in  the  class  of  translation 
invariant  procedures. 

As  mentioned  above,  it  is  difficult  to  imagine  that  one  would  ever  want  to 

consider  adjustment  procedures  outside  the  class  I.   It  should  be  noted 

that  the  class  I  contains  the  class  of  linear  procedures  as  a  proper 

subset;  indeed,  intuitively  speaking,  the  linear  procedures  form  a  "small" 

subset  of  I .  Thus  our  theorem  may  be  described  in  broad  terms  as  dropping 

the  "linear"  restriction  imposed  by  Grubbs  and  Nadler  and  Eilbott. 

Our  application  of  Wolfowitz's  results  can  be  extended  to  a  larger 

class  of  target  destruction  functions,  by  taking  the  loss  function  to  be 

L(y,y)  =  -P  [Procedure  using  S1  =  y(xJ,...,X  ) 
y  n      i      n 

destroys  the  target  on  the  n+lst  round] 

< 

=  -d(y-y(X},...,xl)) 
1      n 

The  X-bar  procedure  will  be  best  (in  I) ,  provided  that  the  destruction 

function  is  a  nonincreasing  function  of   |y  -  y|;  existence  and  finiteness 

of  E_d(-y(xK  .  .  .  .X1))  is  automatic.   (Relationships  between  d(x,y)   and 
u      l      n 

d(x)   are  discussed  below,  in  Section  5.) 
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4.   Summary  and  Conclusions. 

The  foregoing  analyses  indicate  that  an  optimal  (in  several  senses) 
adjustment  scheme  for  the  grenadier  to  pursue  consists  of  shifting  the 
aimpoint  by  an  appropriate  fraction  of  the  miss  location  of  the  preceding 
grenade  burst.   It  is  significant  that  this  procedure  does  not  require 
numerical  evaluation  of  miss  distances,  nor  does  it  require  computations 
with,  or  manipulations  of  the  sequence  of  impact  points  observed  as  firing 
proceeds.  The  procedure  we  have  studied  is  well  known  and  often  used  in 
practice — our  contribution  has  been  to  determine  its  optimality  under 
simplifying  assumptions. 

The  most  stringent  of  the  assumptions  made  in  our  analysis  is  that 
the  grenadier  can  produce  a  shift  in  aim  which  is  precisely  some  fraction 
of  an  observed  miss  distance.  Some  insight  into  the  effects  of  relaxing 
these  assumptions  can  be  gained  by  considering  the  analytical  expression 
in  our  model  for  the  variance  of  aimpoint  after  n  rounds.   It  can  be 
argued  that  the  grenadiers  inability  to  exactly  produce  the  desired  shifts 
in  aimpoint  will  have  the  effect  of  increasing  variance  in  aimpoint. 

If  the  grenadier  makes  normally  distributed  "errors"  with  mean  0 
and  variance  t2,   independent  of  the  actual  miss  distance,  the  overall 
situation  is  as  if  the  ballistic  distribution  has  variance  a2  +  t2  =  1  +  t2, 
with  an  errorless  grenadier.   That  is,  under  the  independence  assumption, 
the  mathematical  model  reduces  to  the  previous  case  with  a  larger  variance. 
The  variance  of  the  aimpoint  for  the  n  +  1 —  round  is  now  (l+t)/n, 
as  compared  to  1/n  with  the  errorless  grenadier. 

A  more  realistic  model  for  grenadier  error  is  that  the  variance 
decreases  with  miss  distance.   Such  would  be  the  case,  for  example,  if 
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the  conditional  distribution  of  grenadier  error,  gives  a  miss  distance 
x,   is  uniformly  distributed  over  the  interval  x  +  px,  where  p  repre- 
sents a  proportion  parameter  measuring  the  grenadier's  skill.  This 
situation  places  a  premium  on  obtaining  rounds  near  the  target,  as  soon 
as  possible.   It  is  anticipated  that  the  X-bar  procedure  will  be  shown 
to  be  optimal  with  a  variety  of  grenadier  error  models  of  this  general 
type;  we  are  continuing  with  work  along  these  lines.   Simulations  with 
models  of  this  type  appear  to  support  a  claim  that  the  variance  in  aim- 
point  does  not  go  down  like  1/n,   as  it  did  in  the  earlier  models  we 

considered.   Rather,  the  variance  appears  to  go  down  somewhat  slower  with 

-3/4 
the  latter  model;  perhaps  as  n      or  so. 

These  variations  in  the  variance  in  aimpoint  under  changes  in  the 
grenadier  error  models  and  parameters  should  provide  a  method  of  assessing 
the  value  of  multiple  shot  grenade  launchers.  This  might  be  accomplished 
by  assessing  the  benefits  of  the  multiple  round  capability  (vice  a  single 
shot  capability)  in  terms  of  the  grenadier  error  characteristics,  which 
in  turn  affects  the  variance  in  aimpoint.  The  aimpoint  distribution  should 
still  be  approximately  normally  distributed,  by  the  central  limit  theorem. 
Hence,  with  a  knowledge  of  the  aimpoint  variance  and  ballistic  variance, 
together  with  a  target  destruction  function,  one  can  compute  the  proba- 
bility of  destruction  with  each  round,  as  well  as  the  expected  number  of 
rounds  required  to  destroy  the  target. 
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5.   Appendix. 

We  have  asserted  that  an  optimal  two-dimensional  procedure  is 
obtained  using  optimal  one-dimensional  procedures  in  each  dimension. 
The  purpose  of  this  appendix  is  to  verify  that  this  is  the  case,  under 
the  following  assumptions: 

(i)   The  range  X  and  deflection  Y  components  of  the  impact  point 
of  a  round  fired  at  aimpoint   (0,0)   are  independent,  symmetric  unimodal 
continuous  random  variables,  with  densities  f(x)   and  g(y),  respectively. 

(ii)   The  target  destruction  function  d(x,y)   symmetric  in  each  argument, 
where  d(x,y)   denotes  the  probability  that  a  round  impacting  at   (x,y) 
destroys  a  target  at   (0,0).   Equivalently,  under  these  symmetry  condi- 
tions, d(x,y)  may  be  viewed  as  the  probability  that  a  round  impacting 
at   (0,0)   destroys  a  target  at   (x,y). 
(iii)   {(x,y):  d(x,y)  >  a}   is  convex,  for  all  a. 

The  approach  we  take  is  to  present  two  "auxiliary"  destruction 
functions,  related  to  d(x,y),  which  are  one  dimensional  and  satisfy 
the  conditions  of  Anderson's  inequality.   It  is  shown  that,  at  each  stage 
of  the  sequential  adjustment  procedure,  optimal  adjustment  for  the 
"marginal"  problems  will  give  the  optimal  "joint"  adjustment.   Suppose 
the  target  is  at   (0,0).   The  probability  of  destroying  the  target,  given 
one  aims  at   (u..  ,u_)   is 

OO    00 

f 

p(u  ,u J  =  P[Target  destroyed|aim  at  (y^y^]  =      d(x,y)f  (x+u1)g(y+y2)dxdy 


00    00 


t 


r 


d(x-p1,y-y2)f (x)g(y)dxdy, 


_00   _00 
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where  f(x+y  )g(y+y  )   is  the  bivariate  distribution  of  impact  points 
with  aimpoint   (y-  ,y?)  .   The  problem  is  to  show  that  this  is  maximized 
with  y   =  y  =  0;   Nadler  and  Eilbott  do  so  [Lemma  2]  using  Anderson's 
inequality,  which  requires  that  d(x,y)   be  symmetric  and   {(x,y):  d(x,y)  >  a} 
be  convex  for  all  a.   Now 


p(ylSy2) 


00   I-  00 


oo   t-  00 


d(x-y1,y-y2)f (x)dx 


g(y)dy 


is  maximized  if  it  is  possible  to  maximize 


d2(y-y2)  = 


d(x-y1>y-y2)f (x)dx 


over  y   for  each  value  of   y  -  y_,   independent  of  the  value  of   y  -  y_, 
and  then  to  maximize 


d2(y-y2)g(y)dy 


over  y  .   Let  d  (x)  =  d(x,y).   As  a  function  of  x  alone,  note  that 
d  (-x)  =  d(-x,y)  =  d(x,y)  =  d  (x) ;   and   {x:d  (x)  >  a}   is  convex,  for 
otherwise   { (x,y) :d(x,y)  >  a}  would  not  be  convex.   Thus  Anderson's 
inequality  implies  that,  for  the  one  dimensional  problem  with  ballistic 
distribution  f(x)   and  any  of  the  destruction  functions  d  (x) ,   the 
maximum  of  d?   is  attained  with  y-  =  0,   independent  of  y_.   Now  apply 
Anderson's  inequality  again  for  the  other  dimension,  with  ballistic  dis- 
tribution g(y)   and  destruction  function  d_(y).   [Note  that 


d2(-y)  = 


d(x,-y)f(x)dx 


00 

f 


d(x,y)f(x)dx  =  d2(y), 
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and  {y:d2(y)  >  a}  =  {y 


d(x,y)f (x)dx  >  a}   is  convex,  again  because 
{ (x,y) :d(x,y)  >  a}   is  convex  for  all  a.]   It  then  follows  that 


d(y-u_)g(y)dy  is  maximized  with  u_  =  0. 

Thus,  given  any  two  dimensional  problem  satisfying  the  conditions 
posed  above,  an  optimal  solution  is  obtained  following  optimal  one 
dimensional  aim  strategies  at  each  stage.   Note  that  radial  symmetry  of 
the  bivariate  destruction  function  d(»,»)   is  not  required. 
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